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ABSTRACT 

A fully Sinc-Galerkin method for recovering the spatially varying stiffness and damp- 
ing parameters in Euler-Bernoulli beam models is presented. The forward problems are 
discretized with a sine basis in both the spatial and temporal domains thus yielding an ap- 
proximate solution which converges exponentially and is valid on the infinite time interval. 
Hence the method avoids the time-stepping which is characteristic of many of the forward 
schemes which are employed in parameter recovery algorithms. Tikhonov regularization is 
used to stabilize the resulting inverse problem, and the L-curve method for determining 
an appropriate value of the regularization parameter is briefly discussed. Numerical exam- 
ples are given which demonstrate the applicability of the method for both individual and 
simultaneous recovery of the material parameters. 


lr rhis research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS 1-18605 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1 Introduction 


In the modeling and control of large flexible structures, one is often required to numerically 
recover one or more material parameters given data measurements at various points. Al- 
though these structures are in general very complex, in many cases the essential features can 
be developed by considering a fixed Euler-Bernoulli beam which is assumed to have Kelvin- 
Voigt damping. In this paper, a fully Sinc-Galerkin method is presented for the numerical 
recovery of the stiffness parameter El and the damping parameter c#/ in the state space 
model 


with 


C(EI } c d 1)u — 0 < x < 1, t> 0 


u(0, <) = u(l, £) = 0, t > 0 

£(»,*) -£<m)- o, <>o 

Sxi 

u(x,0) = -^-(x,0) — 0 < x < 1 

£(EI, c D I)u = -Qp + ( EI ^d^ + CDl Wdx*8t J ’ 


( 1 . 1 ) 


d 3 u \ 


given measurements of tire data at tire points ((xp, ^ 9 ) }p — 1 ... ^ nr tire d o rrrarrr ( 0 , 1 ) x R+. 
From physical considerations, it is reasonable to let the admissible parameter set Q be defined 

by 2 

Q = {(EI,c d I) G n^(°, 1) : EI(x ) >EI o > 0 , c D I(x) > o| 

(see [5]). With this definition, the existence of a unique solution u to the forward problem 
can be obtained on any fixed time interval [ 0 ,t],t > 0 , for / sufficiently smooth. 

The “idealized” parameter recovery problem can then be formulated as follows: determine 
q = (El, cpl) G Q such that 

Cu(-,-,q) = d (1.2) 

where «(•, • ,q ) = C~ 1 (q)f denotes the parameter-dependent state solution and d is the data. 
The observation operator C is given by 


C'P = {ip(x p , t q )} q p= \:: n rZ 


(1.3) 
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and hence Cu represents point evaluations of the solution. Note that (1.2) can be written as 
the operator equation 

)C(q) = d (1.4) 

with the nonlinear operator K given by 

K{q)=CC- l (q)f . 

For reasons similar to those discussed in [11], the problem (1.4) is ill-posed in the sense 
that solutions q (provided they even exist) may not depend continuously on the data d. 
Consequently, some sort of regularization (i.e., stabilization) is required to obtain an accurate 
approximation for q. 

The regularization technique that is used is Tikhonov regularization [13], and the problem 
(1.4) is replaced by the minimization problem 

mmT a (q) (1-5) 

«ieQ 

where 

TM = j{||K(,) - 4“ + cj{q)}. 

Here a > 0 is a regularization parameter which controls the tradeoff between goodness of 
fit to the data and stability. The penalty functional J(q) provides stability and allows the 
inclusion of a priori information about the true parameters. 

Due to the infinite dimensionality both of Q and of the state space, the problem (1.5) is an 
infinite dimensional minimization problem. In order to develop a practical numerical scheme, 
the problem must be replaced by a sequence of finite dimensional problems | that is, one 
must approximate the operator K and minimize the functional T a over a finite dimensional 
admissible subspace of Q. 

The evaluation of JC(q) requires the solution of (1.1). Similar partial differential equa- 
tions must be solved to obtain the components of the derivative fC'(q). The construction 
of an approximate solution to these forward problems commonly begins with a Galerkin 
discretization of the spatial variable with time-dependent coefficients. This yields a system 
of ordinary differential equations which is solved via differencing techniques. For problems 
with nontrivial c D I, it is noted in [l] that the equations are moderately stiff and routines for 
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stiff systems must be employed, thus adding to the expense of the algorithms. This difficulty 
is further augmented by the fact that the time-stepping must be repeated at each step in the 
minimization of (1.5). A final difficulty lies in the need to interpolate at data points which 
do not coincide with the nodes of the ODE solver. 

In contrast, the method of this work implements a Galerkin scheme in time as well as 
space, thus bypassing many of the difficulties associated with time-stepping methods in the 
context of inverse problems. Corresponding results for the heat equation can be found in 
[7], and a detailed discussion of fourth-order results involving the recovery of El in models 
with no damping is given in [11]. 

The fully Sinc-Galerkin method in space and time has many salient features due both 
to the properties of the basis functions and the manner in which the problem is discretized. 
Perhaps the most distinctive feature of the method is the exponential convergence rate 
when solving the corresponding forward problems. Furthermore, the judicious choice of a 
conformal map provides approximate solutions which are valid on the infinite time interval 
rather than only on a truncated time domain. Finally, the discrete system requires no 
numerical integrations to fill either the coefficient matrices or the right-hand side matrix. 
All three features prove to be advantageous when solving the forward problems and hence 
the inverse problem. 

In Section 2, the Sinc-Galerkin system for the forward problem is considered and imple- 
mentation details are discussed. The forward results are then incorporated into the finite 
dimensional minimization problems in Section 3 with the discussion centering around the con- 
struction of the Tikhonov functional. Numerical results are presented in Section 4 along with 
a brief outline of the “L-curve” technique for determining the regularization parameter a. 
Examples are given which demonstrate the recovery of the individual parameters El and 
c D I as well as the simultaneous recovery of both parameters. 
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2 The Forward Problem 


Consider the forward problem 


Cu(x, t) = f(x, t), 0 < x < 1, f>0 

u (0, i) = u(l, t) = 0, t > 0 


du 

dx 


du 
du A 


( 2 . 1 ) 


u(x,0) = — (x,0) = 0, 0 < x < 1 


with 


£«(x, () = ~(x, t) + jL (e/(x)^(x, () + c/(x)^(x, ()) . 

Since a thorough derivation of the Sinc-Galerkin method for problems of this type (with 
cpl = 0) is given in [10], the following discussion contains only that material which is 
needed for the construction of the associated matrix system. 

A fully Sinc-Galerkin method for the approximation of the solution of (2.1) may be briefly 
summarized as follows. For <f>(x) = In (jf- -), T(tu) = ln( w) and positive k x and h t , define 


where 


Si(x) = S(i, h x ) o 4>{x) = sine { ^ X \ 

f*) 

oin/^ 

sinc(x) = 


V h t 


S-(t) = S(j, h t ) o T (t) = sine ( TW~ jh * 

8in(7rx) 


— oo < X < oo . 


( 2 . 2 ) 

(2.3) 

(2.4) 


The basis is then taken to be {5 tJ } with 


= Si(x)Sj(t) } 

and the approximate solution is defined by way of the tensor product expansion 


^ ^ m x = M x -f N x + 1 

U m a m t (2J,0 = 2-/ 


i——M m j~ — M t 


(2.5) 


m t — M t + N t -f 1 . 


4 



The m x ♦ m* unknown coefficients are determined by orthogonalizing the residual with 

respect to the set of sine functions This yields the discrete Galerkin 

system 

(£u m . m ,-/,S p SJ) = 0 (2-6) 

for p = —M x , • • • , N x and q = -M t , • • • , N t . The inner product (•, •) is taken to be 

(F, G) = £ F{x, t)G(x , t)w{x, t)dxdt (2.7) 


where the weight is 

u;(x,t) = u;(xK(0 = (^(x)r J (t(0) _ ^ • ( 2 - 8 ) 


A thorough discussion motivating this choice of weight can be found in [10]. 

The expressions (2.1), (2.7), and (2.8) are then combined to form the system 

Jo J o ^(Um m m t )S p S;WW*dxdt 


(2.9) 


= J°° fS r S',ww'dxit 

for p = — M x , ■ • • , N x and q = —M t , • ■ ■ , N t . 

In anticipation of the parameter identification problems which motivate this analysis, the 
terms El and c D I in (2.9) are expanded as linear combinations of weighted sine functions 
with four Hermite-like algebraic terms added to accommodate the potentially nonzero func- 
tion and derivative values of El and c^I at x = 0 and x = 1. Specifically, this parameter 


basis is taken to be {i>k}k=-M m 


with 


x(l — x) 2 , k — —M x 


M x ) 


(1 — x) 2 (2x + 1), k — —M x + 1 
* V£;(x)Sjt(x), k = —M x + 2, • • ■ , N x — 2 

x 2 (2(l - x) + 1), k = N x — 1 
— (1 — x)x 2 , k = N x . 


( 2 . 10 ) 
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Here Sk(x ) = S(k, h x ) o <p(x) and the basis weight is taken to be 


Ve(x) = w(x) = [x(l — x)]J . 


( 2 . 11 ) 


The finite dimensional approximations for El and c D I are then given by the expansions 

EI m m {x) = Y C k^k{x) ( 2 . 12 ) 

k=-M. 

and 

N m 

co/m.(x) = Y Cki>k(x) . (2.13) 

Jfc=-Af. 

In the forward problem, the coefficients {c fc } and {c*} are known whereas in the corre- 
sponding parameter recovery problems, they are unknown and are determined via methods 
to be discussed in Section 3. The number of basis functions used in the expansions is chosen 
so as to guarantee a square spatial coefficient matrix. 

A quick note should be made concerning the choice of parameter basis and the manner of 
expanding EI mw and cr>I mm . The two derivative-interpolating boundary basis functions are 
added so that these expansions are the same as those used with cantilever or free boundary 
conditions. The choice of (2.11) for the basis weight is certainly sufficient and proves to be 
beneficial when incorporating this forward scheme into a numerical method for solving the 
parameter recovery problem as described in Section 3. 

Sine quadrature is used to evaluate the integrals in (2.9) and hence derive a discrete 
system. For details of the quadrature rule and conditions governing its error bound, see 
[12]; for the purposes of this work, however, it suffices to state the sine quadrature results as 
follows. Let T be (0, 1) or (0, oo) when x = <t> or T, respectively. If F is analytic on T and 
suitably bounded, and if there exist positive constants K , a, and (3 such that 


Hr) 


X\r) 


< K 


c -q|x( t )I j r 6 ^((-00,0)) 


e-'W, reV’UO.oo)) 

where ip = x _1 , then for h > 0 sufficiently small 


(2.14) 


j=-M X(*j) 

The sine gridpoints are given by Zj = ip(jh) — 


< K ie - 2 * d ' h + — e -°Mh + £ e -/SiVfc 
a (3 


(2.15) 
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The expansions ( 2 . 12 ) and (2.13) are substituted into (2.9), and integration by parts is 
used to transfer the derivatives onto the product S p wS*w*. As detailed in [10], the weight 
choice (2.8) guarantees that all boundary terms vanish. The resulting integrals are then 
evaluated via (2.15). For Y — El or Cjjl and 

3>(x) = y(x)-F(0)[(l-x) 2 (2x-f l)]-F(l)[x 2 (2(l-x) + l)] 

— F'(0)[x(l - x) 2 ] - y'(l)[-(l - x)x 2 } , 


the requirement 

|J(x)u(x, t)| < Kx a+ *{ 1 - xf+h^e- 6 ' (2.16) 

guarantees the decay needed to truncate the infinite quadrature rule as specified by (2.14). 
With a ,/?, 7 and 8 specified and M x chosen, the choices 


h x — 



hi — h xy 


a 




a 


M ‘ = l7 Mx + T 


and 

= l^' n + ‘1 

for the stepsizes and summation limits balance the asymptotic errors to at least order 
This rate results from the presence of a sine function in the integral. In 
the above expressions, [•] denotes the greatest integer function. Note that the +1 is unnec- 
essary when or ~M X are integers. 

Given M x , N x , M t , N t and h = h x = h t as defined above, the discrete system for (2.1) is 

AejUCJ + A cd1 UA% + C x UAl = G. (2.17) 


Here 
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and 



where V denotes a diagonal matrix whose entries consist of function evaluations at the sine 
gridpoints. The m x x 77ij matrices U and F are defined componentwise by 


[U]ij = u {j 


and 


[^1‘i — f( x i> tj) • 

It should be noted that the ordering of the coefficients t Hj in U mimics that used in most 
standard time-differencing schemes. This is a matter of convenience since the Sinc-Galerkin 
method is not bound by any specific ordering of the grid. 

As shown in [8], the nit x mi temporal matrices An and An are given by 

a, = + i/] ixmi) , 

/1 “ = [i^ 1 ’ - 1 7 ] 

where I denotes the m* x m* identity, and 1^ and / ^ are defined componentwise by 




0, j = q 


[J«u# = 


IT 

~T’ 

(_ 2)(-l)J'-« 

U - q) 2 

The ra, x m x spatial matrices Agj and A Cd j have the form 


(; - q) 
_2 




3 = q 

» 3 ± q 


^c D I = [$C 2 )X>(p* (J) ) + 2 $< 3 >P(p* (1) ) + S^Vipw)} 
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where 

p*(0 = ¥ ( °c, £ = 0,1,2 , 

= ^ (/) c, £ = 0,1,2 

with c = [c_ Mm ,---,c N .] T and c = [c_ M . , ■ • • , c«J T . Finally, the matrices $ (i) , j = 2,3,4 
and £ = 0,1,2 are defined componentwise by 

[* (i) ]pi = {S p w) {,) (xi) (2.18) 

and 

[¥< q )* = #(xi) (219) 

with the notation on the right-hand sides of (2.18) and (2.19) indicating the j-th and £-th 
derivatives, respectively. The expansions of and in terms of more fundamental 
matrices can be found in [10]. 

For given El and c D I, one then needs to solve the linear system (2.17) for the matrix U. 
As shown on page 414 of [6], (2.17) is algebraically equivalent to 

Au = {C t ® A E i + A tl ® A CdJ + A t2 ® C x } co(U)=co(G ) (2.20) 

where the vector u =co(&) is the concatenation of the m x x m t matrix U obtained by 
successively “stacking” the columns of U, one upon another, to obtain an (m x • m<) x 1 
vector. The system (2.20) can be solved directly via any of the decomposition methods that 
are available for large linear systems. It should be noted however, that for large values of m x 
and m t , the ill-conditioning of the (m x • m t ) x (m x • m t ) matrix A must be considered when 
devising schemes for solving the system (2.20). One facet of current research is directed at 
devising linear algebra techniques which will facilitate the solution of this system. 
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3 The Finite Dimensional Minimization Problem 

As noted in the introduction, the minimization problem 


minT„(,) 

where 

Ul) = Imw - dp + aj(g)} 

is infinite dimensional and thus must be replaced by a sequence of finite dimensional problems 
before a viable numerical scheme can be developed. Following from (2.12) and (2.13), the 
approximating admissible parameter sets are taken to be 

Qm m = {( EI mm ,c D I m . ) ■■ EI mm (x) = Ek'-M. Ck^k(x), 

c D^m m {x) = EkZ-M m CfcV’fc( x )} 

with the basis {ipk} defined in (2.10). For q mm = (Elm*, the associated finite dimen* 

sional optimization problem can then be formulated as 

< si ) 

where 

£,(?».) E 1 {||X(,„.) - dp + aJ(,„.)} . (3.2) 

Note that in solving the minimization problem (3.1), one is actually solving for the vec- 
tors c ~ [c_ w . , ■ ■ • , c Ar»] € , c — [c_m. , • • ■ , cat.] e M m ", or c = [c, c] € iR 2m " depending 

on whether one is solving for EI mm , C£>/ m<) , or both parameters simultaneously. 

With Tip and n q specifying the number of spatial and temporal observation points, respec- 
tively, the approximation K(q mm ) to K{q) is obtained by applying the observation operator C 
in (1.3) to u rnmTnt in (2.5). If the set of observation points {(ip, ^g)}p = i’...’nj can be represented 
as a tensor product of spatial and temporal points, then K(q mm ) has the representation 

^(?mj = C co(U ) (3.3) 

where the matrix U solves (2.17). C is an ( n p n q ) x ( m x -m t ) evaluation matrix which can be 
formulated as follows. Define the n p x m x spatial evaluation matrix E x to have components 

[£x]p,< = Si(x p ), 1 < p < n p , -M x <i<N x 
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and let the n q x m* temporal evaluation matrix Et have the components 

[Et] q .j = S*(t q ), 1 <q<n g , -M t <j< N t . 

Then 

G = Et ® E x . 

It is noted that if the set of observation points is not rectangular as described above, then 
point evaluation can be done directly via (2.5). This latter option is less efficient, however, 
than that defined in (3.3). 

The form of the discrete penalty functional J depends upon whether one is solving for 

EIm w , cplm., or both simultaneously. In the first case, the discrete penalty functional is 
taken to be 

Aim.) = f'[EC m (x)} 7 dx + 5 [\EI mm (x)] 2 dx « c t Qc 

JO Jo 

where the m x x m x matrix Q = Qj + Q j has components 

[Qd]ki « / x)dx, —M x <k,l<N x 

J 0 

and 

[Qf]ki « / M x )M x ) dx , -M x <k,t< N x . 

J o 

The exact representations for the matrices Q d and Q f are given in [11], Similarly, if EI mm is 
known and c D I mtt is unknown, an appropriate penalty functional is 

d (Qm m ) = I [cdI'^Jx)^ dx e f [cD/ mt (i)fii w c T (Jc 

JO Jo 

with Q defined as above. Finally, in the case where both parameters are unknown, the 
discrete penalty functional can be taken to be 

«/(<?m.) ~ c t Qc 

where the 2 m x x 2 m x matrix Q is given by 

Qf + Qj 0 

0 Qf + Q d 
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Although the matrix Q in each case is not sparse, as shown in [11] it is very efficient to 
construct since the component matrices are also needed in the forward solver. Moreover, for 
e > 0 the matrix Q is symmetric and positive definite and hence has a Cholesky decompo- 
sition Q = R T R where R is upper triangular. If c is used to denote c, c, or c, then it follows 
that the penalty term J(q mm ) yields the quadratic form 

c t R t Rc=\\Rc\\ 2 (3.4) 

where || • || denotes the Euclidean norm. This representation for the penalty functional is 
particularly useful both when implementing a scheme to solve the minimization problem and 
when plotting the “L-curve” to determine a suitable regularization parameter a (see [3]). 

To highlight the dependence of the functional T a in (3.2) on the unknown vector c, let 

K(c) = K(q m .) = C co(U{c)) 

where U(c) solves (2.17) for a given c. Noting (3.4), the optimization problem (3.1) can be 
replaced by 

min T a (c) (3.5) 

where 

T.{1) = i{||JT(c) - df + a||flcf }. (3.6) 

To obtain a minimizer for the nonlinear functional T a (c), a quasi-Newton trust region iter- 
ation [2] 

Cfc+l = Cfc + 3 k 

is used. Here s* solves the quadratic programming problem 

min ^ {||/C(c fc ) + K'(c k )s - d|| a + a||72(c fc + s)|| 2 } 

subject to |[a|| < S k with K'(c k ) denoting the Jacobian of K at c k . The trust region radius is 
chosen so that T a (c) has sufficient decrease at each iteration to guarantee convergence to a 
local minimizer of T 0 (for further details about the theory and implementation of the trust 
region algorithm, see [2] or [4]). 
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An important numerical issue in the implementation of the trust region scheme is the 
formulation of the derivative of the operator K . Here the derivative, or Jacobian, is the 
matrix whose I'-th column is given by 

[K\c)\ v = ljm ^[^(c + Te v ) - K{c)\ 
where the standard unit vector e u has components 

r 

1, k = v 

= &vk ~ ^ 

0, k v . 

In the examples of the next section, the Jacobians are calculated with a standard forward 
difference scheme. This scheme is easy to implement and accurate enough for the purposes 
of the method. If further efficiency is desired, a directional derivative scheme such as that 
described in [7] can be used. 


4 Implementation and Numerical Examples 

The three examples in this section demonstrate the use of the Sinc-Galerkin method for 
recovering the individual parameters El and cpl as well as the simultaneous recovery of 
both parameters. In each case the state solution is u(x,t ) = x(l — x) sin(47rx)t 2 e _f and the 
true material parameters are EI(x) = 1 + x + x(l — x) and cpl{x) = 1 -f sin(7rx). For these 
functions, the choices a = (3 = 7 = | and 6 = 1 satisfy the decay condition (2.16). In 
all three cases, the dynamics of the problem are assumed to be modeled by (1.1) with the 
forcing function /(x, t ) consistent with the state solution and the true stiffness and damping 
parameters. 
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In the first example, El is considered known and cpl is numerically recovered using the 
methods of this paper. The roles are reversed in the second example and El is numerically 
recovered while cpl is considered known. The final example demonstrates the numerical 
recovery of both unknown parameters. It should be noted that a comprehensive set of 
examples demonstrating the application of the Sinc-Galerkin method to models of undamped 
beams ( cpl = 0) can be found in [11]. 

A very important practical consideration is the choice of the regularization parameter a 
for a given (error-contaminated) data set. If the error in the data is discrete and random, 
then under certain conditions the method of Generalized Cross Validation (GCV) can be 
used to determine a suitable value of a (see [14]). A second method for determining the 
regularization parameter is to plot the norm of the penalty functional, ||J2c a ||, versus the 
norm of the residual, ||A'(c e ,) — d|| (see [3] or [9]). Here c a denotes the solution to (3.5). In 
this way, one can qualitatively get an idea of the compromise between the minimization of 
these two quantities. The scheme for determining the “optimal” regularization parameter 
consists of finding those values of a such that (||/f(c a ) — rf||, ||i2c Q ||) lies in the “corner” of 
the resulting curve, known as the L-curve. The use of this technique for determining suitable 
choices for the regularization parameter is demonstrated in the examples. 

In each of the three examples, the data was sampled on a regular grid {( x p , t,)} C 
(0, 1) x (0,3]. Nine equally spaced points x p = pAx,Ax = .1, were taken in space and six 
equally spaced temporal points t q = qAt , At = .5, were taken for a total of n = 54 data 
points. To the data, we added a pseudo-random noise vector c from a Gaussian distribution 
with mean 0 and standard deviation a chosen so that the noise-to-signal ratio £r/||d|| = .001; 
that is, noise = .1% of the signal. 

The summation limits were taken to be M x = N x = M t = 8 and N t = 4 as dictated by 
the choice of decay parameters. Hence m x = 17 basis functions were used in the expansion 
of each material parameter. The startup vectors in each example were chosen to reflect the 
positivity and general endpoint behavior of the true parameters. All problems were run with 
sixteen place accuracy on a Vax 8550. 
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Example 4 A. 

In this example, the parameter El is considered to be known and the parameter cpl 
is to be numerically recovered. The startup vector was taken to be the m x x 1 vector 
c 0 = [2, .5, .5, • ■ • , .5, .5, — 2] T . For varying values of the regularization parameter a, the 
L-curve is plotted in Figure 4.1. Note that the points a = 10 5 and a = 5 X 10 6 yield 
points (||if(c a ) — d\\ y ||7?cJ) in the “corner” of the curve. For a = 5 x 10“ 4 , ct = 10 -5 and 
a = 10~ 7 , the plots of the true and approximate damping parameters are given in Figure 4.2. 
It can be seen that the “corner” value, a = 10 -5 , provides a good choice for the regularization 
parameter whereas a = 10“ 7 is not large enough to damp out the contribution due to the 
smaller singular values. Finally, the choice a = 5 x 10 -4 causes too much smoothing and 
information about the parameter is lost. The results from this example demonstrate the 
viability of the method for problems in which El is known and cpl is unknown. 



Residual Norm || /T( cv,) — d\\ 


Figure 4.1: The Tikhonov L-Curve for Example 4.1. 
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Figure 4.2: True and Approximate Damping Parameters for Example 4.1 

(a = 5 x 10~ 4 ), (a = lO' 5 ), •••(« = 10~ 7 ), (True). 


Example 4-S 

Consider now the case where the stiffness parameter El is considered to be unknown 
and the damping parameter cpI is assumed to be known. Here the m x x 1 startup vector 
c 0 = [1, 1, .5, • • • , .5,2, 0] T was used. Since the L-curve for this example was very similar 
to that in Example 4.1, the true and approximate stiffness parameters corresponding to 
a = 5 x 10 -4 , a = 10 -5 and a = 10 -7 were computed with results given in Figure 4.3. It is 
again noted that the “corner” value of a = 10 -5 provides a good choice for the regularization 
parameter whereas the choice a = 5 x 10 -4 causes too much smoothing. Finally, the error 
contributions due to the smaller singular values become more apparent with a = 10 -7 . 
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2.2 



Figure 4.3: True and Approximate Stiffness Parameters for Example 4.2 

(q = 5x lO" 4 ), (a = 10~ 5 ), ■ • ■ (a = 1(T 7 ), (True). 


Example 4-8 

In this example, the method is used to simultaneously recover both the stiffness parameter 
El and the damping parameter cpI . Following from the previous two examples, the 2m x x 1 
initial vector had the form cq = [co, Co] T . The true and approximate material parameters 
corresponding to a = 10~ 5 are plotted in Figures 4.4 and 4.5. From the figures it is noted that 
although the method is accurately recovering the shape of the functions, there is more error 
in the magnitude than in the previous two examples. Hence it appears that a larger number 
of state and parameter basis functions are needed to accurately recover both parameters 
simultaneously, and current efforts are directed at devising linear algebra techniques which 
would facilitate the solution of the larger discrete systems. 
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Figure 4.4: True and Approximate Stiffness Parameters for Example 4.3 
(a = 1(T 8 ), (True). 



Figure 4.5: True and Approximate Damping Parameters for Example 4.3 
(a = 10~ 5 ), (True). 
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